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We develop a nonequilibrium mode-coupling theory for uniformly sheared systems starting from 
microscopic, thermostatted SLLOD equations of motion. Our theory aims at describing stationary- 
state properties including rheological ones of sheared systems, and this is accomplished via two 
steps. Firstly, a set of self-consistent equations is formulated based on the projection-operator 
formalism and on the mode-coupling approach for the transient density correlators which measure 
the correlations between the density fluctuations in the initial equilibrium state and the ones at later 
times after the shearing force is turned on. The transient time-correlation function formalism is then 
used which, combined with the mode-coupling approximation, expresses stationary-state properties 
in terms of the transient density correlators. A detailed comparison of our theory is also presented 
with the related mode-coupling theory which is based on the Smoluchowski equation for Brownian 
particles under stationary shearing. 

PACS numbers: 64.70.pm, 61.20.Lc, 83.50.Ax, 83.60.Fg 



I. INTRODUCTION 

Nonlinear rheological behavior of glassy materials un- 
der stationary shearing has attracted considerable atten- 
tion in recent years since it provides additional insight 
into the physics of glass transition P, [3, [3, 0, H, S 
For such systems driven far from equilibrium, the shear 
rate should be regarded as a relevant control parameter 
rather than as a small perturbation 0]. In this paper, 
we develop a nonequilibrium statistical mechanical the- 
ory for glass-forming systems in which the shear rate as 
well as temperature and density can be handled as ex- 
ternal control parameters. This will be done by extend- 
ing the projection-operator formalism 10] and the mode- 
coupling theory (MCT) |TJ to nonequilibrium systems. 

MCT has been known as the most successful micro- 
scopic theory for the glass transition. Indeed, exten- 
sive tests of the theoretical predictions carried out so far 
against experimental data and computer-simulation re- 
sults suggest that the theory deals properly with some 
essential features of glass-forming systems jT3, [l^. It 
is therefore natural that extensions of MCT have been 
attempted to stationary sheared systems. 

At present there exist two different approaches in 
such nonequilibrium extensions of MCT: one based on 
steady-state fluctuations [13, [lB| , and the other based on 
the transient time-correlation function (TTCF) formal- 
ism m, In the former approach, basic objects are 
the steady-state density correlators defined with fluctu- 
ations around the stationary state. Rheological proper- 
ties like the shear stress are then expressed within the 
mode-coupling approximation in terms of these steady- 
state correlators. With the same spirit as MCT for qui- 
escent systems [llj], the structure factor of the sta- 
tionary state, which now depends on the shear rate as 
well as on the wave "vector" q, enters as input into the 



equations describing the dynamics. At first sight, such 
an approach looks quite reasonable, but in fact it pos- 
sesses a conceptual problem. For example, the following 
exact relation holds between the interaction part of the 
steady-state shear stress cr"* and the steady-state pair- 
correlation function gss(r) 



xy d(j){r) 
r dr 



(1) 



for a uniform shear with velocity along the x-axis and its 
gradient along the y-axis |18| . Here p denotes the average 
number density, and 4){r) the pair-interaction potential. 
Since is related to the Fourier transform of 5ss(r), 
Eq. ^ states that S'q and o-gg* should be handled on an 
equal footing, but this aspect is missing in the steady- 
state- fluctuations approach of Refs. [3, EBl where is 
treated as the input while trgg* is the output. In addition, 
it is assumed in Ref. [l3] that the fluctuation-dissipation 
theorem (FDT) holds also in the nonequilibrium station- 
ary state. The use of such an assumption is unjustified 
since the violations of the FDT have been reported in the 
computer-simulation study of sheared systems [3l- 

On the other hand, no such problems arise in the the- 
ory developed by Fuchs and Cates (FC) [l^ [l3| which 
is based on the TTCF formalism [l^ . Starting from the 
Smoluchowski equation for interacting Brownian parti- 
cles under stationary shearing, the FC theory aims at 
describing steady-state properties via two steps: firstly, 
the MCT equations for the transient density correlators 
- the correlators between the density fluctuations in the 
equilibrium starting state and the ones at later times af- 
ter the shearing force is turned on - arc formulated, and 
then the TTCF formalism is used which, combined with 
the mode-coupling approximation, expresses stationary- 
state properties in terms of these transient correlators. 
In this approach, only equilibrium static structure factor 
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is required as input, whereas the steady-state structure 
factor as well as the shear stress a™* are the output of 
the theory. Thus, the aforementioned conceptual prob- 
lem in Refs. 13, H'l does not apply here. Furthermore, it 
is in principle possible using the FC theory to investigate 
the violations of the FDT, although this issue has not yet 
been addressed. 

It is expected on physical grounds that the microscopic 
dynamics does not matter as far as the long-time glassy 
dynamics is concerned. Indeed, it was argued that the 
equilibrium MCT leads to the same glass transition sce- 
nario for both Newtonian and Brownian microscopic dy- 
namics [ill, [3 HSl , and this was confirmed by computer- 
simulation studies [111, [22|. However, it is not a priori 
obvious whether such an equivalent long-time dynamics 
holds true also for nonequilibrium sheared systems. 

In this paper, we develop a nonequilibrium MCT start- 
ing from the SLLOD equations ^] - Newtonian equa- 
tions of motion under stationary shearing - which have 
been widely adopted in simulation studies of homoge- 
neously sheared systems (see, e.g., Refs. [3, H^])- Our 
theory follows the FC formulation in that the MCT equa- 
tions for the transient density correlators are derived 
first, and then the TTCF formalism is used for describ- 
ing stationary-state properties. However, we found that, 
although it is not difficult to adapt the FC formulation in 
Ref. [ll)[i3 to the SLLOD equations at the formal level, 
the resulting equations are too cumbersome to be useful 
in practice. We therefore developed an alternative formu- 
lation to be presented in the following. It is found that a 
new memory kernel enters into our nonequilibrium MCT 
equations reflecting the non-Hermitian nature of the rele- 
vant Liouville operator, which is absent in the FC theory 
formulated with the Brownian microscopic dynamics. In 
what circumstances this additional memory kernel from 
our theory matters is an open question. We shall elabo- 
rate on this at the end of the paper. 

The paper is organized as follows. In Sec. m we de- 
rive exact microscopic equations and relations for sys- 
tems subjected to stationary shearing. These exact re- 
sults serve a basis for the development of our nonequi- 
librium MCT. We will then derive a set of self-consistent 
equations for the transient density correlators based on 
the projection-operator formalism (Sec. IIII[) and on the 
mode-coupling approach (Sec. HV)) . It is then argued in 
Sec. El how the steady-state properties can be evaluated 
within the mode-coupling approximation based on the 
knowledge of the transient density correlators. The paper 
is summarized in Sec. IVIl where a detailed comparison 
of our theory is also presented with the FC theory. Ap- 
pendix [X] is devoted to a summary of miscellaneous ma- 
terials which are necessary in the main text, and to var- 
ious technical manipulations in the derivations of some 
equations. Appendix iBl describes details of the isotropic 
approximation which is useful in practical applications 
of our theory to systems where anisotropy in the density 
fluctuations is small. 



II. MICROSCOPIC STARTING POINTS 

In this section, we derive exact microscopic equations 
and relations subjected to stationary shearing along with 
thermostat. These exact results serve a basis for devel- 
oping a nonequilibrium MCT for sheared systems to be 
presented in later sections. 

A. SLLOD equations of motion 

We shall consider a system of N atoms of mass m 
in a volume V subjected to stationary shearing charac- 
terized by the shear-rate tensor k. For a simple uni- 
form shear with velocity along the a:-axis and its gradi- 
ent along the y-axis, which we consider throughout this 
paper, the shear-rate tensor is = '^dx^S^y with 7 de- 
noting the strain rate. It is postulated that the applied 
shear induces a homogeneous streaming-velocity profile 
u(r) = K • r at position r, assuming that no spontaneous 
symmetry breaking takes place. Newtonian equations of 
motion describing such a homogeneously sheared system 
are the thermostatted SLLOD equations [Toj . 

Yi = — + K- r.„ (2a) 
m 

= 7i- K- Pi - api. (2b) 

Here and p^ refer to the position and momentum of 
the ith particle, = —dU/dvi with the total interaction 
potential U is the conservative force exerted on the zth 
particle by other particles, and api is the thermostatting 
term which prevents the system from heating up due to 
the work done on it by the shearing force. The momenta 
{Pi}, referred to as the SLLOD momenta, are peculiar 
with respect to the streaming velocity u(ri) = • at 
the particle position r^, and satisfy J2i Pi = 0. 

The thermostatting multiplier a controls the kinetic 
temperature or some other quantity such as the inter- 
nal energy. There exist various types of thermostats - 
stochastic or deterministic, and reversible or irreversible 
ones - considered in the literature. Among them, the 
Gaussian isokinetic thermostat has acquired a respected 
status and a special importance [loj . However, from a 
fundamental point of view, there is no privileged thermo- 
stat, and one should not attribute a fundamental role to 
special assumptions about such models since they sim- 
ply describe various ways to take out energy from the 
system. Indeed, it has been conjectured that different 
thermostats may lead to the same steady-state proper- 
ties, in the usual sense of the macroscopic equivalence 
of equilibrium ensembles [l^l . Although no proof is yet 
available, it is at least reasonable to expect that steady- 
state properties do not significantly depend on the types 
of the thermostats, and this has been tacitly assumed in 
simulation studies where various models have been used 
as practical means to control the temperature. 

In the present work, we shall adopt a constant-a ther- 
mostat in which the multiplier a can be regarded as a 
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"friction" constant. As we will see later, this thermostat 
greatly simplifies the equations to be handled compared, 
e.g., to the corresponding equations under the Gaussian 
isokinetic thermostat whose multiplier ac reads (l9| 



(3) 



How the steady-state temperature can be controlled with 
the constant-a model will be discussed in Sec. IIIGI 



B. The Liouville equation 

For nonequilibrium systems described by the SLLOD 
equations, the form of the Liouville equation commonly 
used for Hamiltonian systems should be properly gen- 
eralized to take into account the effect of phase-space 
compression (l9j . The Liouville equation for the nonequi- 
librium phase-space distribution function f(T,t), where 
r = (r^,p^) stands for a phase-space point, is given by 



dt 



r— +A(r)]/(r,t)^-z£t/(r,i). (4) 



The operator is called the /-Liouvillean, and A(r) 
defined by 



(5) 



is referred to as the phase-space compression factor. For 
the SLLOD equations ^ with constant a, one obtains 



A(r) = 



_d_ 



d 



= -3Na. 



(6) 



The formal solution to the Liouville equation ^ reads 

/(r,t)-expH£tt)/(r,0), (7) 

where exp{~iCH) is called the /-propagator. 

The time evolution of phase variables, which by defi- 
nition do not depend on time explicitly and whose time 
dependence comes solely from that of the phase F, is 
determined by 

^^A{T)^t.-^Ai^)^^/:AiT). (8) 

The operator iC is referred to as the p-Liouvillean. The 
formal solution to this equation can be written in terms 
of the p-propagator exp(z>Ci) as 



A{T,t) = exp{ia)A{T). 



(9) 



Let us summarize here for later use relations between 
/- and p-Liouvilleans and corresponding propagators. It 
follows from Eqs. (g]) and ^ that 



(10) 



One can show that iC and iC^ are adjoint operators, and 
this is why the notation iC'' is used for the /-Liouvillean: 

J dr[ijCA{T)]B{r) ^- J dTA{r)[iC''B{T)]. (11) 

This property can be proved from the integration by 
parts. By a repeated use of this property, the following 
relation for the propagators can be derived: 

J dr [e*^*A(r)] B{T) = J dTA(r) [e-'^''B(r)]. (12) 

If the phase-space compression factor A(r) is identically 
zero, then iC'' = iC holds, and the Liouvillean becomes 
self-adjoint, or Hermitian. In general, this is not the case 
for nonequilibrium systems. 



C. Nonequilibrium distribution function 

Let us consider an equilibrium system of temperature 
T to which a constant shear rate j is applied at time 
t = 0, and thereafter the system evolves according to the 
SLLOD equations The p-Liouvillean is given by 



iC = 



(t < 0), 



iCo + iJ~"y + i^a {t > 0). 



(13a) 



Here, an unperturbed adiabatic or quiescent part (iCq), 
a shear part (iC-y), and a thermostat part (iCa) are re- 
spectively given by 



iC - y^fP' "if 

I m dvi * dpi 

i 

iC^ = Y^^{k ■ Ti) ■ - (k ■ p,) ■ 



dpi - 



(13b) 
(13c) 
(13d) 



Since the phase-space distribution function at t = co- 
incides with the equilibrium one, which we choose to be 
the canonical distribution, 

UT) ^ /(r,0) = ^e-f'"^>(^\ Z = J dTe-^"°^^\ 

(14) 

where (3 = l/k^T with /cb denoting Boltzmann's con- 
stant and Hq = pf /2m + U, a formal solution to the 
Liouville equation (H]) for t > is given by 



/(r,i) = e-^^Vcq(r) 



With the identity 



ds I 



(15) 



(16) 
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whose validity can easily be verified by differentiation 
with respect to t, Eq. pS]) can be expressed as 

/(r, t) = /eq(r) + /* ds e~^^'^{-tC^)f,^iT). (17) 
Jo 

Since iCof^qiT) = 0, we get from Eqs. (P. (fTU)) . and (fT^ 

z/:t/^q(r) = */:^/eq(r) + ^£„/eq(^) ~ 3iVa/cq(r). (is) 

The first term in this expression is given by 

z/:^/eq(r) = /3 V[(/« • rO • F, + (k • p,) • /eq(r) 



13k ■ *T/eq(r) 



^/eq(r)a.,. (19) 



Here cr denotes the stress tensor whose elements are 
(^x^^J2^P^p':/^ + r■Ft] (A,/i = x,y,z), (20) 



and in the final equality of Eq. (|T9)) we used the spe- 
cific form KXf^ = jSxx^f.iy for the shear-rate tensor and 
the symmetry — a^x of the stress tensor (see Ap- 
pendix |XT]). The second term in Eq. is given by 



^'^a/cq(r) 



2a 



/eq(r)x(r), 



(21) 



in terms of the kinetic energy K{T) = X^iPf/^*^- We 
therefore obtain 



»^Vcq(r) 



^ /eq(r)a,,(r) + ^/,q(r)<5i^(r), 



(22) 

where we have introduced the kinetic energy fluctuation 
5K{T) defined as 



5K{T) = K{T) - -Nk^T. 
Substitution of Eq. into Eq. (H?]) then yields 



(23) 



/(r,i) = /,q(r)-^ 



2a 



/*dse-^^'n/eq(r)a,,(r)] 

Jo 

/ dse-*^'^[/eq(r)Ji^(r)]. (24) 
Jo 



This expression for the nonequilibrium phase-space dis- 
tribution function plays a fundamental role in the follow- 
ing. Notice that the last term in this expression vanishes 
if the Gaussian isokinetic thermostat [see Eq. ([3])] is used 
for constraining 5K{T) to zero. In this case, Eq. ((24l) 
reduces to the so-called Kawasaki distribution function 
for sheared systems [l^, [l^ . 



D. Transient time-correlation function formalism 

In contrast to equilibrium quantities, the nonequi- 
librium ensemble average {A{t)) of a phase variable A 
depends explicitly on the time t past since the start 
of shearing. Similarly, the time-correlation function 
{A{t -\- T)B{t)*) depends not only on the time difference 
T but also on t. Using the nonequilibrium phase-space 
distribution function /(r,i), {A{t)) and {A{t + T)B\t)*) 
can be expressed as 

{A{t)) = j dr/(r,o)A(i) = j dr/(r,t)yi(o), (25) 

{A{t + T)B{r)*)= j dTf{T,Q)A{t + T)B{i)* 

dTf{T,t)A{T)B{Q)*. (26) 



The two representations in terms of /(r,0) or f{T,t) 
are equivalent because of the relation Hereafter, 
we shall reserve the notation (• ■ ■ ) for representing the 
averaging over the initial canonical distribution function 

/(r,o) = /eq(r): 



(•••) = / dr/eq(r) 



(27) 



It should be remembered, however, that the dynamics 
inside the brackets (• • • ) is governed by the thermostat- 
ted SLLOD equations, and only averages like (^(0)) and 
{A{0)B{Q)*) coincide with equilibrium quantities. 

Substituting Eq. ^ into Eqs. ^ and ^ and then 
using Eq. (fT2| . one obtains 



{A{t)) {Am 



7 

ksT Jq 



ds {A{s)a„j{0)) 



2a 
kBT Jo 



ds {A{s)SK{0)), 



{Ait + T)B{ty) = {A{T)B{Qr) -^[ds {A{s + T)B{sy a^yiQ)) - ^ |J 



(28) 



ds{A{s + T)B{s)*5K{Qi)). (29) 



r 



The expression (|28p relates the nonequilibrium value of a phase variable A at time t to the integral of tran- 
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sient time-correlation function (TTCF) {A{s)(Jxy{0)) - 
the correlation between cr^y in the initial equilibrium 
state, (Jxij{0), and A at time s after the shearing force is 
turned on - and another integral of TTCF {A{s)SK{0)) 
formed with SK{Q). Equation ((29)) is a generalization of 
this TTCF expression to the time-correlation function. 

The system is said to be in a nonequilibrium steady 
state if the ensemble averages of all phase variables be- 
come time- independent. Let us notice that the long-time 
limit of Eq. (pS)) becomes constant if the system displays 
mixing [19| . This feature can be shown by taking a time 
derivative of Eq. 

(30) 

If the system displays mixing [l9| , then all the long-time 



correlations between phase variables vanish. We there- 
fore obtain for t ^ oo 

j^{A{t)) ^ --l-(Ait)){axym 
In 

-j-fiMt))mm^o, (31) 

since the equilibrium ensemble averages (cFxyi^^)) and 
{6K{0)) are zero [see Eqs. !(Ml and (US])]. This indicates 
that the long-time steady state average of an arbitrary 
phase variable becomes constant, i.e., 

\im{A{t)) = {A),,, (32) 

where the steady-state average, denoted by (• ■ ■ )ss here- 
after, is obtained from the t ^ oo hmit of Eq. (f28|) : 



(A),, = {Am ds {A{s)axy{Q)) -— ds {A{s)SK{0)). (33) 

kbT Jo k^T Jq 

Similarly, the t — > cx) limit of {A{t + T)B{t)*) becomes independent of t since the time derivative of Eq. (P^ . 



f^{A{t + T)B{tY) = -^{A{t + r)B{traxym ~ ^{A{t + r) B (t)* S K (0)) , (34) 



becomes zero for i cxd if the system exhibits mixing. The steady-state time-correlation function defined as 



{A{t)B*),s= \im {A{t + T)B{tr), (35) 

t—>oo 



is then given by 



(A(r)S*),, = {A(T)B{0r) I ds {A{s + t)S(s)V„,(0)) - — / ds {A{s + T)B{sySK{0)). (36) 



The TTCF expressions (1551) and ([55]) . relating the steady-state quantities to the integrals of TTCFs, can be considered 
as the generalized Green-Kubo relations (l9| . 

In deriving the nonequilibrium Zwanzig-Mori-type equation of motion to be presented in Sec. IIIIl it is necessary to 
know how the p-Liouvillcan iC behaves inside the time-correlation function. To this end, we first notice from Eq. ^ 

j^{A{t + T)B{ty) = {[iCA{t + r)] B{tr) + {A{t + r) [iCBit)]*). (37) 



Combined with Eq. (IM|) . this yields the desired result 

{iCA{t + T)B{tr) ^ -{A{t + T)[iCB{t)]*) - ^{A{t + r)i?(t)V,,(0)) - ^{A{t + r)i?(t)*5if (0)). (38) 



For systems exhibiting mixing, there holds for t —> oo 

{[iCAiT]] B*),, = -{A{t) [iCB]*),,, (39) 

i.e., the p-Liouvillean becomes Hermitian in the steady 
state. This is expected since the time-translation sym- 
metry is recovered in the stationary state. 



E. Implication of translational invariance 

Since we are dealing with amorphous systems, the equi- 
librium distribution function /eq(r) is assumed to be 
translationally invariant and isotropic. In this subsec- 
tion, it is shown that the nonequilibrium distribution 
function /(F, t) under shear becomes anisotropic, but re- 
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mains translationally invariant. We then discuss an im- 
plication of this property. Our treatment here foUows the 
one presented in Ref. |17| . 

To this end, we shall consider global translation of all 
particle positions, 

r ^ r' where r ■ = r, + a for all i, (40) 

which amounts to the shift a of the coordinate origin. 
Under this shift, the nonequilibrium distribution f{T,t) 
given in Eq. transforms to 

f{T',t) = /,q(r)--l- /*dse-^*(^'>[/eq(rK,(r)] 

Kb J Jo 



2a 



/*dse-'^'(r')^[/eq(r)<5if(r)].(41) 
Jo 



Here we used 



/eq(r') = /eq(r), a.,y{T') - a,y{T), 5K{T') = 5K{T). 

(42) 

These hold since /cq, axy [see Eq. (|A6|], and 6K de- 
pend on momenta and particle separations only. How 
the /-propagator transforms under F — > F' is discussed 
in Appendix lA 21 with the result [see Eq. (|A17p ] 



d_ 
dvi 



(43) 

Here denotes the transposed matrix of k. Because 
of Eq. dUD, we have P/oq(F) = 0, Va^yiV) = 0, and 
VSkIt) = 0, so that 

e-^^'(^''l/eq(r)a.,(F)] = e-^^'^^^V.<,{'C)a.,y{n, 

(44) 

and a similar equation holds in which axy is replaced by 
5K . Therefore, the nonequilibrium distribution function 
f{T,t) remains translationally invariant: 



/(F',t) = /(F,i). 



(45) 



We next consider how the wave-vector-dependent 
phase variable of the form 



(46) 



transforms under the shift of the coordinate origin. It is 
assumed that xf'^ (F) is a function of momenta and par- 
ticle separations only, so that xf'^{T') = xf''^{F). For 

example, X^'^ — 1 for density fluctuations, xj'^ — /m 
for current density fluctuations to be introduced below, 

and X< = p^p':/m - (1/2) Ej^^ir^r'^/rf,)P.i(rrj) for 
the wave- vector dependent stress tensor [see Eq. (|A4[) ]. 
Using the result 



for the p-propagator which is also derived in Ap- 
pendix |X2] [see Eq. (|A16p ]. one obtains 



Aq(F',i) 



„i£(r)t e'^-'^'^-P* ^X^''(F) e^'^' 

i 

s*('i+i-'**>Mq(F,t), 



(r.+a) 



(48) 



where we used xf-'iT') = xf^{r), Pxf^(r) = 0, and 

ga-/«'^-P tgiq-(ri+a) giq-/t-atgiq-(ri-|-a) 

Since the integral over the phase space must agree for 
either integration variables F or F', there holds 

(Aq(t))= I dF/eq(F)Aq(F,t)- J dF'/cq (F' ) (F' , t ) . 

(49) 

Using /eq(F') = /cq(F) and Eq. (gS]), one obtains 



— „i(q+q-'«t)-a 



"^(^qW). 



(50) 



This means that the nonequilibrium ensemble averages 
of phase variables, including steady-state averages, are 
non- vanishing for zero wave- vector only: 



(^q(t))-^q,0(Aq.0(t)). 



(51) 



Similarly, there must hold for nonequilibrium time- 
correlation functions 

(Aq(t + r)Bk(t)*} = J dF/(F,t)^q(F,r)Bk(F,0)* 
= J dF'/(F',i)Aq(F',T)i?k(F',0)*. (52) 

Using Eqs. (gS]) and (gH]), one flnds 

(Aq(t -I- r)Bk(t)*) = e^(q+q-'^-k)- {A^{t + T)i?k(t)*). 

(53) 

This means that Aq{t + T) is statistically correlated with 
B]i(t)* only if k = q(T) with the advected wave vector 
q(T) = q + q ■ K T during the time r, i.e., 

(Aq(t + r)i?k(t)*) =<5k,qM(Aq(t + T)i?q(,)(t)*). (54) 

Thus, as in equilibrium systems, a time-correlation func- 
tion characterized by a single wave vector can be defined: 

C^^(^ + r,i) EE (Aq(i + r)i?q(.)(t)*). (55) 

For the shear-rate tensor nx^ — jSxxS^y, the explicit 
expression for the advected wave vector reads 



q(r) = q + q • KT = {qx, qy + jTqx,qz 



(56) 



Equivalently, one can introduce a time-correlation of the 
following form 



gi£(r')t ^ gi£(r)t ga-K^-Pt 



(47) 



C^^{t + T,t) EE (Aq(_,)(t + r)i?q(t)*). (57) 
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This also follows from Eq. ([53)l by noting that 

q-(I+/«0 = k ^ q = k-(I+Kt)-i =k-(I-KO> (58) 

since the shear-rate tensor satisfies k ■ k = Q. In this 
paper, we shall mainly use the convention (|55p for time- 
correlation functions, and the convention ([57]) will be 
used only for the discussion in Sec. IIII Al 

Finally, we notice for later use the following relation for 
time-correlation functions involving three phase variables 

= -^k.o {A^{t + T)Bq(,)(t)*i?k=o(t)*), (59) 
which can be derived in the same manner as Eq. (|54[) . 



F. Implication of spatial inversion symmetry 

Let us notice that the SLLOD equations ^ are also 
invariant under spatial inversion T — > — F, and hence, 
the /- and p-Liouvilleans have even parity, i>C^(— F) = 
ijC^r) and iC{-T) = i£{T). Since /oq(F), a^y{T), and 
SK{T) also have even parity, so does the the nonequilib- 
rium distribution function according to Eq. (|24p : 



/(-F,i) = /(F,i). 



(60) 



We next consider how the wave-vector-dependent 
phase variable Aq{T,t) of the form given in Eq. ([^5]) 
transforms under spatial inversion. It is assumed that 
Xf-'ir) satisfies 

xf-i~T) ^pAXf--{T) =pAXf-irr, (61) 

where pA denotes the parity of the variable A. Three 
examples introduced below Eq. (|46)) satisfy these rela- 
tions with pp — +1, pjx ~ —1, and Po-^a" = +1- Then, it 
follows from Eq. ^ and i£{-T) = iC{T) that 

Aq(-F, t) = PAA-q{T, t) = pAAq{T, ty. (62) 

Let us consider an implication of Eqs. ([5D|) and ([S^ 
for the time correlation function Cq ^ {t + r, t) defined in 
Eq. (|55p . Since the integral over the phase space must 
agree for either F or — F, there holds 



J dF/(F,t)Aq(F,r)Bq(,)(F, 



0)' 



= yd(-F)/(-F,<)Aq(-F,T)i?q(,)(-F,0)*. (63) 

Using Eqs. ^Q^j and (HH) and noting that J dT ■ ■ ■ = 
M-F)... [e.g., r^dx.-.. ^ /^-d(-x.)--- = 
dxi ■ ■ ■ under Xi — > —Xi], one finds 

C^''{t + T,t)=pAPBC^''{t + T,t)*. (64) 
In particular, autocorrelation function is real: 

C^^{t + T,t) = C^^{t + T,t)*. (65) 
G. Steady-state properties 



Among various stationary-state properties, we shall 
specifically be interested in this paper in the steady-state 



shear stress, kinetic temperature, and density fluctua- 
tions. Here we summarize the TTCF expressions for 
these quantities. 

The steady-state shear stress shall be defined via 



'{crxy)ss/V. 



(66) 



Since the equilibrium ensemble average of axy is zero, 
(cTxyiQ)) — [see Eq. (IA7[) ]. one obtains from Eq. 
the following TTCF expression for a^s- 



7 



kBTV JO 
2a 



ds {axy{s)a„j{0)) 

ds {a,y{s)SK{0)). 



kBTV JO 

The steady-state temperature shall be defined as 

2 



T 



-{K)s 



(67) 



(68) 



in terms of the kinetic energy. Let us show that Tss is 
connected to CTss via a simple relation. To this end, we 
notice that the rate of the change of the internal energy. 
Ho = K + U , is given by 



Ho 



E 



r Pz_Pi 

L TO 



'^a^y - 2aK, (69) 



where the specific form = 'y6\xSi_iy for the shear-rate 
tensor and Eq. (PO)) for the shear stress have been used 
in the final equality. Since there is no internal-energy 
change in the steady state, i.e., (i?o)ss = 0, one finds 
from Eqs. dMl), dMl), and dMl) that 



7 



(70) 



where p — N/V denotes the average number density. 
Thus, it suffices to know cTss to obtain Tgs- This relation 
can also be used to control Tgs by varying the thermostat- 
ting multiplier a. However, a self-consistent treatment 
is necessary in order to set Tss to a desired value, e.g., 
Tss = T which mimics the Gaussian isokinetic thermo- 
stat, since CTss also depends on a. 

In view of Eqs. ([35)1 and ([55)) . the steady-state corre- 
lator F^{t) of the density fluctuations 

Pq(t)EE^e''l-"-'«-7V<5q,0, (71) 

i 

shall be defined via 

F-(t) ^ ^Ihn l(pq(,s + t)pq(,)(s)*). (72) 

One understands from Eq. that F^{t) is a real func- 
tion of time. From Eq. one obtains the following 
TTCF expression for F^%t) 
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= F^{t) - ds (p^is + t)Pq(t)(s)V,,(0)) - ds + i)Pq(t)(s)*<5K(0)), (73) 

in terms of the transient density correlator defined by 

^qW^^(PqWPqW(0)*), (74) 

and other transient cross correlators formed with (Jxy{0) and 6K{0). As a corollary, one gets for the steady-state 
"static" or equal-time structure factor = F^{t — 0) 

/•OO 2 /'OO 



where 5, = Fq{t = 0) denotes the equilibrium static 
structure factor. While Sg depends only on the wave- 
vector modulus, q = |q|, reflecting the isotropy of the 
initial equilibrium state, the steady-state structure factor 
S'q depends also on the direction of the wave vector due 
to the anisotropy of the sheared stationary state. It is 
also clear from Eq. (|75)l that should be considered 

I 



as a dynamic object in the sense it is given by the time 
integrals of the transient time-correlation functions. 

Finally, let us show the connection between the TTCF 
expressions for the steady-state shear stress CTss and the 
structure factor 5^ given by Eqs. (p7|) and ([75]) . respec- 
tively. For this purpose, we first write cTxy{s) appearing 
in the integrand of Eq. (p7|) as [see Eq. (|A6|) ] 



= L y y * T'^ ^'^ V (2^ ^ []^Pq(^)Pq(^) - Ij : (76) 

where we have used f{vij) = J dr f{r)6{r — Vij) and i5(r — r^) = (l/27r)'^ / dqe^"'*-'""'''^-' with r^j — Vj in the 
final equality. Substituting this into Eq. |(S7|1. one obtains 



Here, the definition (|33p for the steady-state average has 
been exploited for the first term, and Eq. ([75)) for the 
steady-state structure factor S'q** has been used for the 
second term. Notice that {S^ — Sg) in the second term 
can be replaced, e.g., by [S^ — l], since isotropic terms do 
not survive after spatial integral involving the anisotropic 
term xy. Equation ([77]) simply expresses that anisotropic 
density fluctuations are responsible for the steady-state 
shear stress. Since the steady-state pair correlation func- 
tion, gss(r), is related to S'q via 

P[9ss{r) - 1] - / e"'''-(^q - 1), (78) 

one understands that the interaction part of Eq. ([77]) is 
equivalent to Eq. ([T|) . Such an equal handling of tiss and 
based on the TTCF formalism is expected since no 



approximation has yet been introduced. 

In the following sections, we will first derive a set of 
self-consistent equations for the transient density cor- 
relators Fq{t) using the projection-operator formalism 
fSec. lIII|) and the mode-couphng approach fSec. lIVp . We 
will then argue that the mentioned TTCF expressions for 
the steady-state properties can be evaluated within the 
mode-coupling approximation based on the knowledge of 
Fq{t) (Sec.|V|. In this way, we construct the nonequilib- 
rium MCT for stationary sheared systems. 



III. ZWANZIG-MORI-TYPE EQUATIONS 

In this section, we derive exact Zwanzig-Mori-type 
equations of motion for the transient density correlator 
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" -iq-K-rst iCt 
dt 



i£s{t)e'''i'^''^'e'^\ (81) 



Fq(t) for a system that is initially at equilibrium and time-deperirfeni, p-Liouvillean defined via 

subsequently subjected to stationary shearing along with 
thermostat. A "standard" approach for a quiescent 
system is that a Zwanzig-Mori equation for a correla- 
tor {A{t)A{0)*) of some phase variable A evolving with 
time- independent p-Liouvillean is derived based on the 
static projection operator onto the subspace spanned by 
A. As we will see below, due to the presence of the time- 
dependent wave- vector advection q(t), this standard ap- 
proach should be appropriately generalized for sheared 
systems. We start our discussion by pointing this out. 



A. Difficulties in applying previous formulations 

Recently, McPhie et al. [1^ developed a projection- 
operator formalism which generalizes the standard ap- 
proach to nonequilibrium systems and allows one to de- 
rive a Zwanzig-Mori-type equation for a transient corre- 
lator {A{t)A{{))*) . Their formalism is also based on the 
time-independent p-Liouvillean and on the static projec- 
tion operator. Using the convention (|57p . one can intro- 
duce the transient density correlator of the form 



(79) 



Thus, apparently, there seems no problem to apply the 
formalism developed in Ref. by setting A{t) = 

Pq(-t){t)- However, pq(_()(t) is not a phase variable 
since its time-evolution is also affected by the wave- 
vector advection, q(— i), and its equation of motion can- 
not be written solely in terms of the time-independent p- 
Liouvillean as in Eq. ([5]). Their formalism, therefore, can- 
not be directly apphed to derive the equation for Fq{t). 
Nevertheless, we mention here that our equations of mo- 
tion derived below resemble those presented in Ref. [1^ 
in that new memory kernels enter in addition to the one 
familiar in the equilibrium Zwanzig-Mori equations. 

More recently, Fuchs and Gates [13] derived the 
Zwanzig-Mori-type equation for -F'q(t), starting from the 
Smoluchowski equation for interacting Brownian parti- 
cles under stationary shearing. It is not difficult, at least 
formally, to adapt their formulation to the SLLOD equa- 
tions, and we briefly summarize here its consequences. 

Because of the equivalence of the particles, the tran- 
sient density correlator Fq(t) defined in Eq. (fM)) can be 
written as 



1 



^q(0 = ^(PqW(0)Vq(i)) 



(80) 

where = e''^ ''^ denotes the density of a single tagged 
particle (labeled s), which is identical to the others. Here- 
after, the absence of the argument t implies that the 
associated quantities are evaluated at t — 0. By this 
trick of singling out a particle, the motion of the collec- 
tive density fluctuations pq can be described by one, but 



Based on the p-Liouvillean iC for the SLLOD equations, 
the operator iCs{t) can be worked out explicitly, 

iCs{t) = iC — iq ■ K ■ r s + iq ■ K ■ (ps/m)t. (82) 

Integrating Eq. ((8T|) in time, one obtains 

^qW = (pr-?'"'^^^Vq)- (83) 



Here e+ denotes the time-ordered exponential, where ear- 
lier times appear on the right. This expression also ex- 
plains why the formalism developed in Ref. (26j , which is 
based on the tiiae- independent p-J^ionviWean, cannot deal 
with Fq(<). 

Equation (|83p can be handled by manipulations based 
on the static projection operator Vs = Pq)(l/S'q)(pq*, 
and one can derive the following exact Zwanzig-Mori- 
type equation of motion for Fq(t) in the same manner as 
detailed in Ref. 'XL]- 



d 



d_ 
dq 



q-K- — Sq 



Fq{t)+ dsKqit,s)Fq{s) ^0. 
Jo 

(84a) 

Here the memory kernel is given by 
1 



+ Qs«'Cs(i )Pq ),(84b) 



X e 



in which Qg = I — Vs with / denoting the identity op- 
erator. Equations ([M]) serve as the starting equations 
for Brownian particles exhibiting overdamped dynamics, 
since in this case the velocity entering into Qsi£s{t)Pq is 
proportional to the force, and hence, K(t, t') essentially 
describes the fluctuating- force correlations. Such an in- 
corporation of the fluctuating-force correlations is essen- 
tial in developing self-consistent equations for Fq{t). 

On the other hand, we need an additional Zwanzig- 
Mori-type equation for K{t,t') in constructing a self- 
consistent theory since the time derivative of the (pe- 
culiar) momentum is proportional to the force in the 
SLLOD equations For this purpose, one needs to 
introduce a tmie- dependent projection operator onto the 
subspace spanned by QsiCs{t)pq. We found that the 
resulting equation of motion for K{t,t') is too cum- 
bersome to be adopted as our starting equation. (See, 
e.g., Ref. [131 for the application of the time-dependent 
projection-operator formalism.) 

Here we shall take an alternative route. Adopting the 
original definition 



1 



1 



^qW = ^(PqWPq(t)(0)*) = ^ ( [e'"Pq] ) , (85) 



N 
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of the transient density correlator, we will first derive 
an exact continuity equation which relates Fq{t) to the 
transient cross correlator 



1 



1 



= ]^OqWPq(t)(0)*) = ^([e'^^J^] P*^^) , (86) 

for A = x,y,z, formed with the current density fluctua- 
tions jq defined by 



J_„jqri 



(87) 



Here it is necessary to take into account all the A compo- 
nents of jq due to the anisotropic nature of the sheared 
system. We will then derive a Zwanzig-Mori-type equa- 
tion of motion for H^{t), which can be done via a partial 
use of the static projection operator as we will see below. 



B. Continuity equation 

We start with the time-evolution equation for the num- 
ber density fluctuations. Since [see Eqs. fT^ ] 



d 



(88) 



one finds the following continuity equation for the 
sheared system relating the partial time derivative of 
Pq(i) = e»^Vq toj^(t) = e*^*j^: 



d_ 
di 



q • K 



9q 



Pq(i) = 'iq-jq(0- 



(89) 



Likewise the partial time derivative of the density fluctu- 
ation at the advected wave vector, p*^^^Y is given by 



d 



d 



dt ^ 



Pq(t) = 0' 



(90) 



in deriving which we have noticed that the shear-rate 
tensor satisfies k k = 0. One can readily obtain from the 
above two equations that the transient density correlator 
^q(^) ^ (Pq(OPq(t))/^ the transient cross correlator 
H^{t) = {j^(.t)P*ciit))/N is connected via 



d 



d 



dt ^ 



Fq(t) =^q•Hq(^). 



(91) 



One therefore gets for [t] = e'"^* jq 



d_ 

dt 



q • K • 



9q 



j^(t)=e'^*z/:oJ^-[/fjq(t)]^-aj^(t). 

(93) 

It is straightforward to obtain from this equation and 
Eq. {nOD for the correlator H^{t) ^ P* u-.) / N 



d 



d 



dt ^ dc\ 



^q^''^f"q(i)/ 



~[/^.Hq(<)]^-ai/^(t). (94) 



It is already clear at this point that one cannot de- 
rive a closed equation for the "longitudinal" component 
q • Hq(i) alone due to the presence of the second term 
on the right-hand side of Eq. ([94]) . This reflects the 
anisotropic nature of the sheared system. We also notice 
that the thermostatting multiplier a can be taken outside 
of the ensemble average in the last term of Eq. ((94|) since 
we have adopted the constant-a model. If, for exam- 
ple, the Gaussian isokinetic multiplier q;g were used [see 
Eq. ([3])], then one would have to consider an additional 
equation of motion for (l/^)([e'^*aGjq] Pq(t))- Thus, a 
considerable simplification is achieved via the adoption 
of the constant-a model. 



D. Projection-operator formalism 

In the following, we shall apply a projection-operator 
formalism, but only to the first term on the right-hand 
side of Eq. (|M)) . As will be shown below, this can be 
done via a static projection operator, and thereby the 
aforementioned difficulty connected with Eq. (|84bp can 
be avoided. In this way, we complete the derivation of the 
Zwanzig-Mori-type equation of motion for (t) , which 
together with the continuity equation (j9ip provides our 
starting equations for developing a nonequilibrium MCT 
for transient density correlators. 

For this purpose, let us introduce the static projec- 
tion operator V onto the subspace spanned by pk and 
j-;^ (At = x,y,^). Since (p kP^) ^ 5k,k'iV5fe, (j^j^^;) = 
'5k,k''^A^i^w^ with V = \/k^Tjm denoting the thermal 
velocity, and {pkjj^*) — (remember that the averaging 
is over the initial canonical distribution), the projection 
operator V is given by 

^^ = E(^^^)^Pk+EE(^^T)^Jk^ (95) 



C. Exact equation for the transient cross correlator 



We next derive an exact equation of motion for the 
transient cross correlator H^{t). We start from 



i£jq = (i£o + + ^-Cq) 



d 



iCoj^ +q K- —j^ - {k ■ jq)^ - aj^. (92) 



The complementary projection operator is defined by 
Q = I — v. One can easily show that V and Q are 
idempotent and Hermitian. 

The time evolution of iCpj ^ ap pearing in the first term 
on the right-hand side of Eq. (|94p shall then be separated 
into parts correlated and uncorrelated with {pk,j^}- 



iCt 



„iCt 



jct 



Q^-CoJq. 



(96) 
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As derived in Appendix I A 3[ one obtains 



On 



(97) 



In the right-hand side of this equation, we have used 
^QiCtQ „ giQCQtQ Tff]^[(.]^ holds due to the idempotency 
of the operator Q. Let us introduce 



and hence, the first term on the right-hand side of 
Eq. lIM)) is given by 



with 



Kit) 



R = QiCfij = iCfij - iqx—Pq, 



(101) 



(102) 



For the second term on the right-hand side of Eq. ([M)) . 
we use the identity 



where we have used Eq. ^7^. Since = Qjq(*) = 0, 

there holds 



JCt ^ ^QiCt ^ ^iC(.t-s)-p^J^^Q^Cs (gg) (i?^(Op;(t)) = and {Rl{t)f^;-) = 0. 

^0 



(103) 



to obtain 



Thus, i?q(i) is always uncorrelated with {pkjj^}, and we 
follow the usual convention to call this phase variable the 
random or fluctuating force. 

In terms of the fluctuating force Rq{t), the second term 



/ dse'^^^~'''>ViCe'^^^'QiCQj^. (100) in Eq. (UnOl) can be expressed as 
Jo 



£ ds e^('-^^V^CR\{s) = £ ds Y,{[^CRl{s)] + £ ds J] ^([z/:i?^(.)] ^■^*)^e'A*-)j-^ 



k A" 



iC{t-s) -M 



q(s)' 



(104) 



where the last equality holds since {R\{s)f^) is nonzero only for k = q(s) [see Eq. ((54)) ]. We also noticed that the 
ensemble averaged terms are independent of the phase and are unaffected by the propagator. The evaluation of the 
ensemble averaged terms in the integrands of Eq. p04p is presented in Appendix lA 41 and the results are given by 



([z/:i?^(s)]p;(,)) = 



kBT 



Ms)' 

2a 



(105) 
(106) 



Let us notice that Eq. (|106p has been simplified due to the adoption of the constant-a model for the thermostat [see 
the comment below Eq. (|A32p ]: otherwise, e.g., when the Gaussian isokinetic thermostat is used, one has to add a 
term (i?^(s) S[aGJq(*)]) to the right-hand side of Eq. piM|) . 

With Eqs. UnOll-dMl) and pin)) -pi} ^ . we now obtain 



JCt 



7^£dsi:,^^(.)e^^(*-^)j^(,) +a£ds*iVq^(s)e^^(*-'^)pq(,) -a^^^ (107) 
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Here we have introduced the following memory kernels: 



1 



m 



'q(t)/ 

A/ 



2 



9m 



(108) 
(109) 
(110) 
(111) 
(112) 



Substituting Eqs. (HH), ([Ml), and (fTOT)) along with Eq. (fT03)) into the first term on the right-hand side of Eq. we 
finally obtain the following Zwanzig-Mori-type equation for H^{t): 



d 



d 



dt dq 



H^it) = iqx^F^it) - [k ■ Hq(t)]^ - aH^it) - J] / ds M^f^is) H^^^^{t - s) 



7 ds zL^is) (t - s) - 7 ^ ^' L;^'^(s) H^^^^ {t - s) 

+ aj\stN^is)F^^,^it^s)~aY,J^dsN^^'^{s) H^^^^ (t-s). (113) 



Here, we have noticed (1/A^)([e*^(* ''Vq(s)] Pq(i)) — 
Fq(,)(i-s) and (l/A^)([e^^(*-^)j^(,)] p*^^,^) = H^^^^{t-s). 
One can easily confirm that these are consistent with the 
definitions ^ and ((86l) . 

The memory kernel M^^{t) describing the fluctuat- 
ing force correlations is already familiar from the equi- 
librium Zwanzig-Mori equation of motion for the density 
correlator 11]. The additional memory kernels L^{t) and 
L'^^^{t) are due to coupHngs between the fluctuating force 
and the shear stress, and N^{t) and N^'^^lt) are associ- 
ated with couplings between the fluctuating force and 
temperature fluctuations. In the following section, we 
introduce mode-coupling approximations for these mem- 
ory kernels to obtain a set of self-consistent equations of 
motion for the transient density correlators. 



IV. MODE-COUPLING APPROXIMATION 

We have encountered five memory kernels in the 
Zwanzig-Mori-type exact equations of motion for the 
transient correlators. We need to invoke approxima- 
tions for these memory kernels in order to obtain closed 
equations for Fq{t). In this section, we apply the mode- 
coupling approximations [llj | to these memory kernels. 

The basic idea behind MCT is that the fluctuation of 
a given dynamical variable decays, at intermediate and 
long times, predominantly into pairs of hydrodynamic 
modes associated with quasi-conserved dynamical vari- 



r 



ables. It is therefore reasonable to expect that the decay 
of the memory kernels at intermediate and long times 
is dominated by those mode correlations which have the 
longest relaxation times. The sluggishness of the struc- 
tural relaxation processes in glass-forming systems sug- 
gests that the slow decay of the memory kernels is basi- 
cally due to couplings to pair-density modes. The sim- 
plest way to extract such a slowly decaying part is to 
introduce another projection operator V2 which projects 
any variable onto the subspace spanned by PkPp, i.e.. 



r.x = j:ixpip;)^p,p,. 

k>p 



(114) 



Here we already used the static version of the factoriza- 
tion approximation introduced below [see Eq. (|120p ]. It 
is readily verified that 'P2 is idempotent and Hermitian. 

The first approximation in the mode-coupling ap- 
proach thus corresponds to replacing the propagator 
giQCQt governing the time-evolution of the memory ker- 
nels by its projection on the subspace spanned by the 
pair-density modes, e'^'''^* « 7'2e*^'^^*7^2- Under this 
approximation, the memory kernel M^^{t) defined in 
Eq. (fT08)) is given by 



1 



1 



me 



iQCQt 



;[e'2^2*P2<]P2i?^(V. (115) 



The expression for the projected random force 7^2^^^ is 
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derived in Appendix I A 51 within the convolution approx- 
imation for triple correlations, 

(PqPkPp) ~ S^M+pNSqSkSp, (116) 

and is given by 




k>p 



Here Cq is the direct correlation function defined via 

pCq = l^^. (118) 
Jq 

Substituting Eq. (fTTTll into Eq. pTSjl . we obtain 

M^''{t) = ^J^Yl '5q,k+p^q(t),k'+p'[fcACfc -I-PACJ 

k>pk'>p' 

X [k'^ck' +p'^Cp,W^'-^'p^Pp]pt,,p;,). (119) 

The final approximation in the mode-coupling approach 
is to factorize averages of products, evolving in time with 
the propagator e^^''^*, into products of averages formed 
with the variables evolving with e^^* (factorization ap- 
proximation) : 

= 5k',k(t) V,P(t)^'^k(i)^pW- (120) 

Here the translational invariance of the sheared system 
is taken into account [see Eq. ([54| ]. Applying this ap- 



Here S'^ = dSg/dq. It is anticipated that this memory 
kernel becomes relevant only if significant anisotropy is 
developed in the density fluctuations. This is because the 
shear stress axy entering into its defining equation p09p . 
which is reflected in the quantities in the second square 
brackets in Eq. (I122|) , is intrinsically an anisotropic quan- 
tity. For example, one finds from Eq. ()122|1 that L^{0) = 
reflecting the isotropy of the initial equilibrium state. 

The other memory kernels defined in Eqs. (|110p - (|112|) 
are found to vanish under the mode-coupling approxima- 
tion as demonstrated in Appendixes lA 71 and lA 81 

L'^^^it) ^ 0. (123) 
N^{t) = 0, Afq^^(t) = 0. (124) 

Thus, only the memory kernels M^^{t) and L^{t) survive 
under the mode-coupling approximation formulated with 



proximation to Eq. (|119p . we obtain 

2 p 

M^'^it) - J dk[kxck+pxcp] 

X [/c^(t)cfe(t) +P^(t)cp(t)]Fk(t)Fp(t), (121) 

where the wave vector p in this and the following expres- 
sions for the memory kernels abbreviates p = q — k, and 
should not be confused with the momentum variable. 

In the absence of shear, the MCT expression (|12ip re- 
duces to the one familiar from the equilibrium MCT [llj 
describing nonlinear interactions of density fluctuations, 
called the cage effect, relevant for structural slowing 
down. The matrix structure as well as the wave-vector 
dependence in M^^{t) can be simplified, i.e., it can be de- 
composed into longitudinal and transversal components 
which depend on the modulus q = |q| only, and this 
is possible because of the isotropic nature of the quies- 
cent equilibrium system. In the presence of shear, on 
the other hand, the "dephasing" of the vertex function 
in Eq. p2ip occurs, which reduces the nonlinear interac- 
tions, and hence, enhances the structural relaxation. In 
addition, the structure of Mq^(t) cannot be simplified in 
a mentioned way due to the wave-vector dependence of 
the vertex function and of the transient density correla- 
tors, which are associated with the anisotropic nature of 
the sheared system. 

The memory kernel L^{t) defined in Eq. (fTog]) can be 
handled in a similar manner under the mode-coupling 
approximation, and its detailed derivation is presented 
in Appendix I A 61 with the result 



Fi,{t)Fp{t). (122) 

I 

the projection operator 7^2 ■ 



V. STEADY-STATE PROPERTIES 

In this section, we provide the TTCF expressions for 
the steady-state properties (see Sec. IIIG|) under the 
mode-coupling approximation. This enables one to ob- 
tain the stationary-state properties based on the knowl- 
edge of the transient density correlators Fq{t). 



A. Remarks on TTCF expressions 

Let us first notice that the transient time-correlation 
functions appearing in the TTCF expressions in Sec. Ill Gl 



kxkyjt) _^ PxPyjt) '^p(t) 

Ht) Skit) pit) Spit) 
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can be abbreviated as 

Gx{t) = ([e'"X] a^y), Hx(t) = ([e^"X] 5K). (125) 

For example, the TTCF formed with ^^^^^(O) in Eq. ^ 
is given by {[e'^'~^crxy\<^xy) : and the one in Eq. ([751) by 

As discussed in Appendix IA 91 the functions Gx{t) and 
Hx{t) evolve in time within the subspace orthogonal to 
{Pk,ik}' i-^-' there hold 



Gx{t) = {[i 



iQCQt 



QX] Qaxy), 



Hx{t) - {[e'^^^'QX]Q6K), 



(126a) 
(126b) 



in terms of the projection operator Q complementary 
to V defined in Eq. ((95)) . This feature is exactly the 
one sheared with the memory kernels [see Eqs. (jlOip . 
(fT03l) . and (fT08l) - (fTT2l) ]. Thus, no additional approxima- 
tion than those introduced in Sec. lIVl is necessary to deal 
with Gx{t) and Hx{t). The only difference here is that, 
since both axy and SK are "zero wave- vector" quantities, 
the second projection operator V2 given in Eq. pi4p has 
to be replaced by V2 defined via 



k>0 



(127) 



Thus, under the mode-coupling approximation, only 
those contributions abbreviated as Gx{t) survive in the 
TTCF expressions for the steady-state properties. 



B. Steady-state shear stress 

With the results in the previous subsection, the TTCF 
expression (|67p for the steady-state shear stress under 
the mode-coupling approximation is given by 



^^^^fc^y ds([e^e^e.^o^^^]pO^^^). (131) 
Substituting Eq. (|129p into this expression yields 



y ^k' 



k>Ok'>0 

x([e^2£S.^^^*]^^,^*,). (132) 



Applying the factorization approximation (|120p . one gets 
([e'^^^VkPyPk'P^) 



[e^^Vklp^ 



[e'^-'pl] Pk' 



4',k(.)A^'j^k(s)2 



(133) 



We thus obtain under the first mode-coupling approxi- 
mation, in which the propagator e^^^^* is approximated 
by the projected one V^e'^^^*^V^, 

Gx(t) = {[e'^^^'vlQX]V°^axy), (128a) 
Hx{t) = {[e'^'^^'V°2QX]V^5K). (128b) 

Here we have noticed Qdxy — c^y and Q5K ~ SK [see 
Eq. (ITTS)) ]. 

The evaluation oi'P2crxy is presented in Appendix IA 101 
with the result 



where in the final equality we have noticed that F\^{s) 
is a real function of time [see Eq. ([S5)) ]. This leads to 
the following MCT expression for the steady-state shear 
stress CTss in terms of the transient density correlators: 



2(2^ 



^Wrfk%Mf)||(f)^,(,)^ 

KK{S) Jk'^k(s) 



(134) 

The steady-state kinetic temperature Tss can then be ob- 
tained via Eq. (|70)) . 



'P2'^xy 



kj^T ^ ^ kxky Sf^ 

k>0 



N ^ k SlP'^'"'^- 



(129) 



k 



In view of Eq. ([231) . one easily understands that 
{5Kp\^pl^ = 0, and hence. 



(130) 



C. Steady-state density fluctuations 



With the remarks in Sec. IV Al the TTCF expression 
([75)) for the steady-state density correlator F^{t) under 
the mode-coupling approximation is given by 



7 



Nk^T 



ds{[e^^''^^Q{p^{t)p*,,,}]V^,axy) 



(135) 



Here we do not apply V2 to Q{Pq(i)Pq(f)} since it already has the form of the density product. Let us notice that, since 
Pk^Q — [see Eq. d7T|) ] and j^^q = (1/m) X^iPf = Ne the comment below Eq. (pbl) ]. it follows from Eq. (1551 ) that 



(/'q(*)Pq(t) Pi) = '5k,0 (Pq(i)Pq(t) Pk=o) 



and (Pq(<)p*(,)jr) 



^k,o (Pq(i)/3q(4) jk=o) " 0. Ouc therefore obtains 
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'^{/Oq(*)/Oq(f)} = [see Eq. ([95])], and hence, Q{Pq(OPq(t)} = Pq(*)Pq(t)- Thus, we have for the integrand of Eq. (|135p 

loo k 

where we have used Eq. (|129p . Here we apply the factorization approximation [see Eq. (|120p ] 

([e'^2£S.^^(i)^* ^^^] ) ^ ([e^^Vq(0] Pt) ([e'^^V^fi,] Pk) = iV'^^q(t + s)Fq(,) (s), (137) 

where in the final equality we have noticed that -F'q(f)(s) is a real function of time [see Eq. ([55)1 ]. This yields the 
following MCT expression for the steady-state density correlator F^{t) in terms of the transient density correlators: 

F-{t) = Fq(0 + 7 r ds^-^^f^^F^it + .)Fq(,)(.). (138) 

Jo q(t + s) 

As a corollary, we obtain for the steady-state structure factor 5'^'^ — F^{t = 0) 

S-^S, + j ds'-^^F,{sr. (139) 

I 



Let us see the connection between Css and un- 
der the mode-coupling approximation. By comparing 
Eqs. and one finds 

where we have noticed that the isotropic term in 5^*^ does 
not contribute to the integral. Thus, ags and are 
handled on an equal footing naturally under the mode- 
coupling approximation. Compared to Eq. (|77p . the ki- 
netic part is missing here since only the interaction part 
is dealt with under the mode-coupling approach. In addi- 
tion, since SyS^ = pc'f. [see Eq. ()118p ]. the bare potential 
in Eq. ([77)1 is replaced by the "rcnormalized" [1Cp| direct 
correlation function in Eq. (|140[) . Finally, we notice that 
Eq. (I140|) can directly be derived from Eq. (|129|) by ap- 
proximating (Tss K, —{V2<Jxy)ssly and using the defini- 
tion5^^ = (l/7V)(pkPj;)ss. 

VI. SUMMARY AND DISCUSSION 

In this paper, we developed a nonequilibrium MCT 
for uniformly sheared systems starting from microscopic, 
thermostatted SLLOD equations of motion. Our the- 
ory aims at describing stationary-state properties includ- 
ing rheological ones, and this is accomplished via two 
steps. Firstly, a set of self-consistent equations of motion 



is formulated for the transient density correlators Fq{t) 
based on the projection-operator formalism and on the 
mode-coupling approach, which enables the calculation of 
Fqit) provided the static structure factor Sq of the initial 
equilibrium state is given as input. The transient time- 
correlation function formalism is then used which, com- 
bined with the mode-coupling approximation, expresses 
stationary-state properties in terms of iVj(t). Thereby, 
steady-state quantities like the shear stress (Tss, temper- 
ature Tss, density correlators F^{t)^ and the structure 
factor can all be calculated in terms of Sq. We also 
addressed how the steady-state temperature Tgg can be 
controlled using the constant-a model for the thermostat: 
this can be done via a self-consistent treatment of the 
thermostatting multiplier a based on Eq. ([70)) . Our the- 
ory is able to treat CTss and on an equal footing, which 
is missing in the steady-state-fluctuations approach of 
Refs. [l3, EBj- In addition, we need not assume the va- 
lidity of the fluctuation-dissipation theorem in sheared 
states, which was necessary in Ref. [l3 |. 

The transient density correlators Fq(<) thus play a dis- 
tinguished role in our approach. Let us collect here all 
the relevant MCT equations for Fqit) derived in Secs. lIIII 
and llVl to highlight new features of our theory compared 
to the equilibrium MCT (Tl| and to the nonequilibrium 
MCT developed by Fuchs and Cates (FC) for sheared 
Brownian systems [1^ [l3] ■ The exact Zwanzig- Mori- type 
equations consist of the continuity equation 



d_ 
dt 



q ■ K 



d_ 



Fq(t) =Zq.Hq(i), 



(141a) 
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and the time-evolution equation for the transient density-current cross correlator H^{t) 



dt ^ dq 



- ^ f dsM^'^{s)H^, {t-s) + j f dstL'^is)F^^,^it-s). (141b) 
„ Jo Jo 

In this equation, we already omitted those memory kernels which vanish under the mode-coupling approximation (see 
Sec. IIVI) . The MCT expressions for the memory kernels M^^{t) and L^{t) are given by 



pv 



2(27r)= 



dklkxCk + PxCp][kf_,{t)ck{t) + Ptj.{t)cp[t)]F]^{t)Fp{t), 



(142a) 



2(27r)3 



dk [ kxck + pxcp ] 



kxky{t) S'm^ PxPy{t) ^p{t) 



m 

I — 



Pit) 



p(t) 



Fkit)Fpit). 



(142b) 



Here p = q — k. Compared to the equilibrium MCT [Tl| , 
new features entering here are (i) the replacement oid/dt 
by [d/dt — q • /t • (d/dq)], (ii) the presence of the second 
(shear) and the third (thermostat) terms on the right- 
hand side of Eq. (jl41bp . (iii) the matrix structure of 
the memory kernel AI^^{t) describing the fluctuating- 
force correlations which cannot be decomposed into the 
longitudinal and transversal parts, and (iv) the pres- 
ence of the additional memory kernel Lq(i). Further- 
more, when compared with the FC theory |l6l . IT^ . we 
see in addition to those rather trivial differences reflect- 
ing the Newtonian and Brownian short-time microscopic 
dynamics (v) the memory kernel in the FC theory de- 
scribing the fluctuating-force correlations, to be denoted 
as Afq*-^(t, t'), has a different mathematical structure in 
that it depends on two times t and t' after the shear- 
ing force is turned on, while only one time enters into 
our M^^(t), and (vi) the memory kernel corresponding 
to Lq(t) is absent also in the FC theory. 

The first and second features just mentioned arise from 
the shear part (iCj) and the thermostat part (iCa) in the 
p-Liouvillean for the SLLOD equations [see Eqs. (fT5|) ]. 
which are absent in the p-Liouvillean for quiescent sys- 
tems. The third feature reflects the anisotropic nature of 
the sheared system: in the presence of shear, the longitu- 
dinal and transversal current density fluctuations cannot 
be separately handled as in isotropic systems since their 
cross correlators do not vanish. In this connection, we 
notice that the second term on the right-hand side of 
Eq. (|141b[) . which cannot be expressed in terms of the 
A-component H^{t) alone, also reflects the anisotropy 
of the sheared system. Therefore, without introducing 
any further approximation (see below), Eqs. p41ap and 
(|141bp cannot be combined to yield a single second-order 
integro-differential equation for Fq(t) as in the equilib- 
rium MCT [ll|. The fourth feature originates from 
the non-Hermitian nature of the p-Liouvillean describing 
nonequilibrium dynamics [see Eq. (|38p ]. i.e., the pres- 



ence of the additional memory kernel L^{t) is expected 
on general grounds. 

The fifth feature, when compared with the FC the- 
ory, is due to different strategies employed in deriving 
the Zwanzig-Mori-type equations for Fq{t): the two-time 
structure in Mq'~^{t,t') is an exact consequence of the 
Zwanzig-Mori-type equations ([51)) for Fq{t) upon which 
the FC theory is based (see Ref. [l3l)i while the one- 
time structure in our M^^{t) follows from another exact 
equation (|94p to which the projection-operator formal- 
ism is applied (see Sec. IIII D[) . One therefore cannot 
judge which of the memory kernels is superior at the 
formal level: we can only state that ours has a simpler 
mathematical structure. Furthermore, both the memory 
kernels Mq^{t) and Mq'~^{t,t') under the mode-coupling 
approximation describe essentially the same physics con- 
cerning the competition between the cage effect and the 
shear advection of density fluctuations (see also below). 

Thus, the sixth feature mentioned above, i.e., the pres- 
ence/absence of the memory kernel Lq{t) is the major 
difference between our and the FC theory. It is not likely 
that this difference originates from the different micro- 
scopic dynamics - Newtonian or Brownian - adopted in 
these theories since, as we stated above, the presence of 
such a memory kernel is expected on general grounds. 

It is anticipated that the memory kernel Lq(i) becomes 
relevant only if significant anisotropy is developed in the 
density fluctuations. This is because the shear stress (Txy 
entering into the defining equation (|109p of L^it) is in- 
trinsically an anisotropic quantity. We indeed confirmed 
from our preliminary numerical calculations based on the 
MCT expression (jl42bp that the contribution from Lq(i) 
is quite small under the isotropic approximation for the 
density fluctuations to be discussed below. It will be in- 
teresting to pursue in what circumstances this additional 
memory kernel becomes important whose presence is nat- 
urally expected for nonequilibrium sheared systems. 

To further facilitate the comparison of our theory 
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with the equihbrium MCT and with the FC theory, 
the MCT equations pH]) and pl^ shall be simpli- 
fied using the isotropic approximation introduced in Ap- 
pendix IB II Such a simplifying approximation is also 
useful in practical applications of our theory to systems 
where anisotropy in the density fluctuations is small. 

The MCT equations ((BT3)) . (|BI9)) . and (|B23|) derived 
in Appendix IB II under the isotropic approximation shall 
be rewritten in the following form for the normalized 
transient density correlators (pqit) = Fq{t)/Sq: 



relaxator equation 



7^^ / ds/^^°(s)05(,)(t 



0. (143) 



Here all the functions depend on the wave-vector mod- 
ulus only; the dot denotes the partial time derivative; 
fi^ = q'^v'^/Sq the square of the characteristic frequency 
relevant for the short-time dynamics; and q{s) = q[l + 
(7s)^/3]^/^ the modulus of the advected wave vector un- 
der the isotropic approximation. The memory kernels, 
from which is factored out following the convention 
in the equilibrium MCT [ll[, are given by 



mf (t) = 
Tit) = 



(7) 



q,k,pW0fcW0p(i), (144a) 
{t)Mt)<Pp(t), (144b) 



with the time-dependent vertex functions 



[q-kck+q- pcp] 



V 



X [q-kcfc(t) + q-pcp(t)], (144c) 
(7)' _ SqSkS„ 



q,k,p 



3 ^1 + (7<)V3 2(27r)-V 



^ [q ■ kcfe + q • pCp] 



C/ CI 



^k(t 



it) 



s, 



(144d) 



The resemblance of these equations to those in the equi- 
librium MCT [ll[ is apparent: the major differences are 
the dephasing in the vertex function V^^p(t) for 'm'^°{t), 
which enhances the relaxation of the density fluctuations, 
and the presence of the additional memory kernel 1^^° (t) . 

Now, let us "derive" the MCT equations for sheared 
Brownian systems, starting from Eq. p43p with the pro- 
cedure adopted in Ref. ^ for converting the microscopic 
dynamics from Newtonian to Brownian. Assuming that 
the "friction" constant a is large, the inertia term in 
Eq. p43p shall be neglected. As a result, the general- 
ized oscillator equation (|143p is specialized to generalized 



(l)qit) + Tq<j>q{t) + Tq / ds m'^'° (s) ^q^s){t " s) 



+ 7r, / ds /;^°(s) [t - s) = 0, (145) 
Ja 

where we have defined Tq = ^q/a. This equation, com- 
bined with Eqs. p44ap and (I144cp and neglecting l'-^°{t) 
which is found to be small from our preliminary calcu- 
lations, is formally identical to the corresponding equa- 
tion in the FC theory. [See Eqs. (4)-(6) of the second 
article cited in Ref. [l^. There is a minor difference 
that <t>q{s) {t ~ s) at the advected wave number q{s) enters 
into the third term in Eq. (|145p . while 4'q{t — s) at the 
wave number q appears in the corresponding FC equa- 
tion. Again, this reflects the difference of the starting 
Zwanzig-Mori-type equations.] In this sense, our and the 
FC theory are equivalent. But, it should be remembered 
that this holds only under the isotropic approximation: 
when anisotropy in the density fluctuations is significant, 
one has to go back to Eqs. (|14ip and (|142p . and the pres- 
ence/absence of the memory kernel L^{t) may have sig- 
nificant consequences. 

Finally, we notice that the anisotropic nature of the 
steady-state quantities can still be captured within the 
mentioned isotropic approximation for the transient den- 
sity correlators. This issue is discussed in Appendix IB 21 
In particular, it is argued there why the steady-state 
shear stress CTss can be evaluated under the isotropic 
approximation, as done in the application of the FC 
theory [IQ], although this quantity should vanish under 
isotropic density fluctuations. 
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APPENDIX A: MISCELLANEOUS MATERIALS 
AND DETAILS OF SOME DERIVATIONS 

This appendix is devoted to a summary of miscella- 
neous materials which are necessary in the main text, 
and to various technical manipulations in the derivations 
of some equations. In these derivations, we repeatedly 
use the relation 



{{zCoA)B*) = -(A{iCoBy), 



(Al) 



which holds for the unperturbed or quiescent p- 
Liouvillean iCo given in Eq. (fT3b|) . and 



d 



dA 



(A2) 
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where F.^ = —dU/drf denotes the A component of the 
conservative force acting on the ith particle. These re- 
lations, well-known from equilibrium statistical mechan- 
ics [To'], hold here since the averaging (• • • ) in this paper 
is defined with the canonical distribution function [see 
Eq. ()27p ]. Also, terms involving odd number of momen- 
tum variables vanish under such canonical averaging. 



1. Microscopic expression for stress tensor 

Here we summarize the microscopic expression for the 
stress tensor. For simplicity, we deal with quiescent equi- 
librium system for which the p-Liouvillean is given by 
iC = iCo [see Eq. (|13bp ]. In handling sheared systems, 
momenta appearing in the following expressions should 
be understood as pecuhar or SLLOD momenta [l9| . 

The wave- vector dependent stress tensor tjq^ is intro- 
duced via the continuity equation for the current density 
fiuctuation = J2iiPi/''^) exp(iq • r^) 



^CoJ^^Y. — <'^ (A3) 



and is given by [T3| 



Here r,- 



1 „M 

-L \ ^ ij ij 



T.lp'p'/^-^T.^P^i'^^^) exp(zq.r,). 



1 — exp(— iq • r) 



rv', and 



Pq(r) =r0'(r)- 



zq • r 



(A4) 



(A5) 



in which 4>{r) denotes the pair-interaction potential. Ob- 
viously, cTq'^ is a symmetric tensor. 

The "stress tensor" referred to in the main text is the 
zero- wave- vector limit of fi^^: 



Exploiting the isotropy of the quiescent equilibrium sys- 
tem, one can show that [lOl 



(A7) 



The equivalence of the expression (|A6|) and Eq. ([^ 
in the main text can be demonstrated as follows. Since 



J 



Fi = X^j^^i -f^u where F^j denotes the force acting on the 
ith particle from the jth particle and Fj; — — F^^ due to 
Newton's third law, there holds 

i i j^i j i^j 

= ^EE^^^F^^-- 

i j^i 

Expressing the force Fy in terms of the pair-interaction 
potential as 



UL i I ij 



(A9) 



one obtains 



(AlO) 



indicating the equivalence of Eqs. (|A6|) and (^0)1 . 

2. Propagators under the global translation 

Here we discuss how the /- and p-propagators trans- 
form under the global translation F ^ F' defined by 
Eq. (1^0]) . The p-Liouvillean corresponding to the SLLOD 
equations is given by [see Eqs. p3)) ] 



r d 
iC{T) = y^^^{pi/m+K-ri)- — + {Fi-K-pi-ap^) 

i * 

and it transforms under F ^ F' to 



_d_ 
dp 

(All) 



iC{T') = iC(T) + a • k'^ • P with P = E 7^ ' ^^^'^^ 



since Pi and F^ = J2j^i ^ij (where F^j denotes the force 
acting on the ith particle by the jth particle and is a 
function of only) are not affected by F — > F'. Here 
denotes the transpose of k. 
Let us notice that, when iC{T) acts on a phase variable 
Ar(F) that depends on momenta {pi} and particle sep- 
arations {vij} only, there holds PX(T) = 0. Therefore, 
the only term in iC(T) that does not commute with P is 
the second term in Eq. IjAlip . for which we have 



P.- 



(A13) 
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We therefore obtain 



P.z£(r) - ^CiT)P. = [Mn • r.)] ■ |- = E (E ^^.^0] |x = E E 5^' 



i A 



and hence, there holds 



(a • • P) i/:(r) - i£(r) (a . • P) = Y.axKj^[P^iC(T) ~iCir)P^] 

X,fi 

= E^^'^A^^EE'^^'A'^ 

= E <^A • '«)a'A P\' = 0, 

A,A' 



(A14) 



(A15) 



since the shear-rate tensor satisfies k k — 0. Thus, ijC{r) 
and a • k"^ • P commute. This means that /-LiouviUean 
i£^(r) and a • ■ P also commute since the difference 
between /- and p-Liouvilleans for the SLLOD equations 
with the constant-ck model for the thermostat is simply 
a constant [see Eqs. ^ and PH]) ]. 

Using the Campbell-Baker-Hausdorff theorem which 
states that e-^'^'^ = e-^ for commuting operators A 
and B, one obtains from Eqs. (jA12|) and (jAlSp 



JC{T)t+a-K^ -Pt _ iC(T)t a-K^ -Pt 



Similarly, there holds for the /-propagator 



(A16) 



relation iCopq = iq • jq, the former is given by 

= -^q.kO'^ (*q • j;)) = <5q,k tqxNv^ (A18) 



For the latter, we use Eq. (|A3[) to obtain 



since only odd number of momentum variables are in- 
volved. It thus follows from these results and Eq. ((95|) 



(A17) 'p^CoJ^^J2(^^CoJ^]pt)J^^p^^^qx^Pq. (A20) 



3. Derivation of Eq. (fOT]) 



Here we derive an expression for ViCoj^. To this end, 
one needs to evaluate the ensemble averages {[iCojq] /oj.) 
and {[iCojq]ji^*) [see Eq. ([95])] • Using Eq. dH]) and the 



4. Derivation of Eqs. ([105]) and (fT06l) 



Using Eq. ((38)) . the ensemble averages in the integrands 
of Eq. (|104p are given by 



Since Qi?q(s) = Rq{s) [see Eqs. (|10ip and (|102p ] and the operator Q is idempotent and Hermitian, the above equations 
can be written as 

{[zCR^qis)] p^j,)) - [QiCpq^s)]*) - J^iKi^) QKis)'^^y]) - QKis)m), (A23) 

mKi^)] & = -iKi') [Q^^fqi.P - ^(K(') QKis^'^^y]) - ^(Ki^) Ql&K])- (A24) 
In the following, we will show that 

(R^is) [QiCpqis)]*) - 0, (A25) 

[Q^CJ^^,^V) = {Kis)R'^l,^). (A26) 
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Substituting these results into Eqs. (|A23|) and (|A24p yields Eqs. (jlOSp and (|106p . respectively. 
To derive Eq. (|A25p . we first notice from Eq. 



i£pq(5) = iq{s) • jq(,) +^i[q(s) • n ■ r^] e^'^^'*^''"^' = iq(s) • jq(,) + q- k 



d_ 



(A27) 



since k ■ k = 0. We therefore obtain, since Qj^, 



d 



On the other hand, it follows by taking a partial time derivative of the first relation in Eq. (|103p that 



(A28) 



d_ 



d_ 



JQCQsj^X 



(A29) 



where in the final equality Eq. (j90p has been used for the second term. The first term in this expression is zero since 
Q is Hermitian and QPq(^s) ~ ^-'^^^ obtain 



0. 



(A30) 



Equation (|A25P then f ollow s from Eqs. (|A28P and (IA30p . 

We next derive Eq. (|A26p . To this end, we notice from Eq. ([9^ that, since k • k = 0, 



d 



q(s)- 



We therefore obtain, since Qjq(s) = 0, 



q(s)J 



Q ■ K ■ —- 7 , X 



(A31) 



(A32) 



Notice that the thermostat term ctjq(s) does not contribute here since we have adopted the constant-a model. If, e.g. 
the Gaussian isokinetic thermostat is used, the contribution (i?q(s) Q[Q!Gjq(*)]) cannot be discarded. 

The vanishing of the second term in Eq. (IA32P can be demonstrated as follows. Using the following equation 



d 

^jq(^i - q • • — J, 



d_ 



(A33) 



a partial time derivative of the second relation in Eq. (|103p is given by 



{{ 



ds 



}<(:)) + 



= ([Q^/:i?^(s)]J^*)) + (i?^(^ 
I 



q ■ K 



• d 
d 



dq 



Jq(s) 



(A34) 



The first term is zero since Qjq(s) — 0, and this leads to 



(A35) 



One therefore obtains from Eqs. (lA32l) and (jASSl) 

- {<i^)Kls))^ (A36) 

where in the final equality we have used Eq. (|102p for 
the definition of the fluctuating force. This completes 
the derivation of Eq. (|A26p . 
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5. Derivation of Eq. (fTTfl) 

Here we calculate the projected random force P2Rq- 

V2RI = Y.Kp*^pI)n^P^Pp- 

k>p ^ 

To this end, we need to evaluate [see Eq. (|102p ] 

{R^pIpD = ([*^oJ^] PkPp - ^qx^^{P^PtPl)- (A38) 

Oq 

Using Eq. (|A1[) and the relation iCoPq ~ iq-jq, the first 
term is given by 

([i^oiq] pIp*p) = -Oq [i^oPk]* P*p) - (jq pl [i^oPp]* ) 
= {jqi^-3t)pl) + {jqPli^P-rp)) 

iNv^[kxSp+pxSk]. (A39) 

For the second term in Eq. (jA38[) . we use the convolution 
approximation pi6p : 

^2 

iQ\-E-{p^PtPp) ~ Sq.k+piNv'^ qxSkSp. (A40) 

Oq 

One thus obtains from Eqs. (|A38|) - (|A40|) 

(RqPtPp) = -Sq^i^+piNpv'^SkSp[kxCk +P\Cp], (A41) 

in terms of the direct correlation function [see Eq. piSp ]. 
Substituting this result into Eq. (jA37| finally yields 




k>p 

I 



6. Derivation of Eq. (fT22|) 

Here we derive the MCT expression for the mem- 
ory kernel Lq(t) defined in Eq. (|109p . Under the first 
mode-coupling approximation e*^-^^* « 7'2e*^'^^*'P2 (see 
Sec. IIVP , one obtains 

(A43) 

Since 7^2 -Rq is already given in Eq. (|A42p . we only need 
to consider 



'P2Q[Pq{t)(^xy] = ^({Q[Pq(t)ga:y]}PkPp) ^25 5 P^PP' 
k>p ^ P 

(A44) 

Let us start from Q[Pq{t)0'xy], for which we need to 
know the averages {[Pq(t)0-xy]pt) and {[pq(t)0-xy]j^*) [see 
Eq. The latter is zero, {[Pq(t)0'xy]ji^*} = 0, since 

this term involves odd number of momentum variables 
only. For the former, one obtains using Eq. (pO| 

{[Pq{t)'^xy] Pk) = '5k,q(t)(Pq(t)CTa;yP*(t)) 

since the kinetic-part contribution from axy vanishes. 
Applying Eq. (IA2p to this equation yields 



d 



([Pq(t)^..]Pk) = -'5k,q(t)A:Br(5^x,— e'^i^*)-'---') 



'5k,q(t) Nk^Tqy[t)j^^ = - Jk.q(t) Nk^T'i^^S'^^,), (A46) 



since qx{t) — qx [see Eq. and dSq(^t)/dqx = (fe/'z)'S'^(t) ■ We therefore obtain from Eq. ([55)1 



and hence 



Now let us calculate 



:P[Pq(t)Cr:jv] = y'([Pq(t)CTxi;] Pk)T77rPk = -feB^^^^^ --2^Pq(t) , (A47) 



Q[Pq(t)^xy\ - Pq(O^X. + fcBT^^^Pq(,). (A48) 

q Oq(t) 



mP^it)^xy]] Ptp*p) = {[Pqlt)CTxy] plp*p) + ^B^^^ ^ (Pq(t) pJIpP • (A49) 
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Using Eq. ((20l) . the first term is given by 



i j e m 

d 



= -'5q(t),k+pA:BT( x,^(e''i(*)--'e-^'''-^e-^P-^-)), (A50) 

i,j,i,m 



where we have employed Eq. (jA2p in the second equality. The calculation of this term can be continued in the same 
manner as in Eq. (|A46[) with the result 

{[PqW'^xy] PIP*p) = -^q(t),k+P^Br(fcy^ +Py^)(Pq(t)P^pP. (A51) 

Substituting this into Eq. (IA49P yields 

{[Q{P^it)a.y}]plP;) = -'5q(t).k+pfcBT{fc, A +p^A. _ MM^hp^^^.^p*^p*)_ (A52) 
This expression can further be simplified under the convolution approximation (|116p . 

(Pq(t)PkPp) ~ Sqit)M+pNSq(^t)SkSp. (A53) 

Let us notice that, when q(t) = k + p, there holds 



Similarly, we have 



-g;^[Sq{t)SkSp] — -^^S'q^^^-^SkSp + ^Sq(t)SkS'p. (A55) 



It then follows from Eqs. (|A52p - (|A55|) that 



([QW)-..}] pIpD = -5^(t)M+pNk^TSq^,)SkSp[^^^ + ^|}, (A56) 



and substituting this into Eq. (IA44|) yields 

P2Q[p,Waxy] = -^^,(0 E ^.w..+p{^| + }p.cPp. (A57) 

k>p ^ 



Substituting Eqs. (|A42p and (|A57|) into Eq. (IA43[) and then using the factorization approximation (|120p . we finally 
obtain with p = q — k 

kxky{t) S'f.f^t) PxPy{t) ^p(t) 



L'^{t) = --—, d\.[k,ck+p.cp] F^(t)Fp{t). (A58) 

I 

7. Derivation of Eq. (|123p mation formulated with V2 ■ We start from 

(A59) 

under the first mode-coupling approximation e^^^Q* « 
Here we show that the memory kernel L'^^{t) defined 7'2e*^'^^'7'2 (see Sec. In the following, we demon- 

in Eq. (fTTU|) vanishes under the mode-couphng approxi- strate {{Q[3'^(^t)'^ xy]} p'iP*^) = 0, i.e., p2Qbq(t)CTj;y] = 0, 
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which completes the derivation of L'^'^{t) = 0. Eq. ((95)) ]. The former is zero, {[j^(^t)0-xy]pt) = Oj since 

Let us start from Q[j!^(^f-^<^xy], for which we need to this term involves odd number of momentum variables 
know the averages {[j^(^t)Cr^y]p^) and {[j!^i^tfxy]jk*) [see only. Using Eq. the latter reads 



In this equation, the kinetic-part contribution survives 
only when (i) i — j = £, fj, = v = y, and (ii) i = 
j = £, fj, = y, I' — X, and the potential-term contribution 
survives only when i = £, fi ^ v. This leads to 

since the potential-term contribution vanishes after ap- 
plying Eq. (|A2|) . We therefore obtain from Eq. ([M)) 

k 

= mv^[S^xjl^f^+S^.yf^^t)i (A62) 

and hence 

Since the right-hand side of this equation involves odd 
number of momentum variables only, there holds 

{{Q\J^it)^^y]}plP;) = Q- (A64) 

8. Derivation of Eq. (fT24l) 

Here we show that the memory kernels N^{t) and 
N^^>'{t) defined in Eqs. (fTTTjl and fm]) vanish under the 
mode-coupling approximation formulated with 7^2- We 
start from the following expressions under the approxi- 
mation e*2£St p^f,iQCQtj)^ (ggg Sec.HVl: 

Kit) « ^^^^l^{[e^Q^Q^V2Rl{t)] V2Q[p*^(t)6K])^ 

(A65) 

(A66) 

In the following, we demonstrate ({Q[/Oq(t)^i^]} P^Pp) — 
and {{Q[j'^i^t)SK]} pIpI) = 0. This means 
V2Q.[p^(t)5K] = and r2Q[j^(^f^SK] = 0, and hence, 

completes the derivation of N^{t) = and N^^t' it) = 0. 

Let us start from Q[pq(t)<5/'ir], for which we need to 
know the averages {[Pci{t)^K]p'^) and {[Pq{t)^K]j^*) [see 
Eq. ([M]) ]. In view of Eq. ([^5]) . one easily obtains 

([Pq(t)'5i^]Pk) - 0, ([Pq(*)5i^]jT) = 0. (A67) 



I 

Thus, V[pq(^t)^K] — 0, and hence, 

Q[p^it)SK] = p^it)SK. (A68) 
Likewise, one obtains 

QKwm=f^^t)SK. (A69) 
It is then obvious that 

{{Q[p^it)SK]} p*^p;) = 0, {{Q[j^^,^SK]} pip;) = o. 

(A70) 



9. Derivation of Eqs. (fT26l) 

Here we show that the functions Gx (t) and Hx (t) de- 
fined in Eq. (|125p evolve in time within the subspace 
orthogonal to {/Ok, j^}, i.e., there hold 

Gx{t) = ([e*^*^]fT,;,) = ([e'^'^^^Q^lQfT^^;), (A71) 
Hxit) = ([e*^*X]5i^) = {[e:'^^'^'QX]QSK), (A72) 

in terms of the projection operator Q complementary to 
V defined in Eq. (j95p . Before embarking on the deriva- 
tion, let us notice 

QcTxy = cr^y and Q6K = SK. (A73) 

The first relation follows from (cTxyPk) = 
S]i,o{(yxyPk=o) = [see Eq. ([Til)] and {(Jxyii^*) = 0, and 
the second relation can be derived in a similar manner. 
Thus, the presence of the operator Q in front of a^y and 
6K in Eqs. (|A7ip and (|A72p is irrelevant. 

In the following, we shall deal with the function Gx{t) 
only, since Hx{t) can be handled in a similar manner. 
Applying the identity 

e'^' = e'^St _^ f dse*^(*-") i£V e'^^\ (A74) 
Jo 

[notice the difference in the order of operators compared 
to the identity ([M)) ]. one finds 
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Jo 

= e'^St^ + T.jk;l ds{[^^^^X]pl)e^^('-^hCp^ + E E ]^ ds{[e^^^'^X]jry'^'-'^^Cj^, (A75) 

I 



where we have used the definition ([95]) of the operator V 
and noticed that the ensemble averaged terms are inde- 
pendent of the phase and are unaffected by the Liouvil- 
lean and the propagator. Let us notice here that 

iCpk and iCjj^ ^0 for k ^ 0. (A76) 

The former is obvious in view of Eq. ([88]) . while the lat- 
ter can be derived on the basis of Eq. ([5^ by noticing 

tCoj^ = (l/m)!:^*^^^^ [see Eq. dMI] and j^^o = 
{^/'m)J2iPi = [see the comment below Eq. ([2bl) ]. 
Equation (IA76P simply expresses the fact that the den- 
sity and the current density, the latter being defined for 
sheared systems in terms of the peculiar momenta, are 
conserved variables. 

Now let us consider the transient correlator 
{[e^^*X]axy) formed with the "zero wave-vector" 
quantity axy The translational invariance implies that 

([e^^(*-^Uk] <Jxy) = <5k,o ([e'^(*-^^Ak=o] a.^,) = 0, 

(A77) 

for — iCpi^ and iCj^ because of Eq. (|A76p . Thus, 
there is no contribution to ([e*^*A'] axij) from the second 
and third terms on the right-hand side of Eq. (|A75P , and 
we obtain 

Gx{t) = {[e'^'X] axy) - {[e'^^'X] a^y). (A78) 
Since the operator Q is idempotent and Hermitian, one 



finds using Eq. (|A73p 

Gx{t) = ([Qe'^2*X] Qa.,) 

= {[Qe'^''^'QX]Qaxy) 

= {W^''^'QX]Qnxy), (A79) 

where in the second equality we have noticed 

Qe^^S* = Qe'SCQtg^ (A80) 

This completes the derivation of Eq. (IA7ip . and 
Eq. (|A72p can be derived in a similar manner. 

10. Derivation of Eq. (fT29l) 

Here we derive the expression for P2(^xy For this pur- 
pose, we need to know the average {uxyp-kPi^)- Using 
Eq. (pni) . this average can be written as 

= -fcBrE(^^^M), (A81) 

where we have used Eq. (jA2p . Since dpu/dyi — ikye^^'^^ 
and dp\^/dkx = i XiC*'' ''', we obtain 



(CT^yPkPk) = -kBT{iky(^{^Xi(^'^''^^pl^ -iky(^{^x,e '''■'■')pk)} 

i i 

I 



It then follows from Eq. p?7)) that 



V^Gxy ^ ^{oxy P\,pl)j^j^p\^pt 

k>0 

= -—l^^^P^P^- (A83) 

k>0 ^ 



APPENDIX B: ISOTROPIC APPROXIMATION 



In this appendix, we shall introduce the isotropic 
approximation which considerably simplifies the wave- 
vector-dependent MCT equations (|14ip and (|142p for the 
transient density correlators. Such a simplifying approx- 
imation is useful in practical applications of our theory 
to systems where anisotropy in the density fluctuations 
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is small. We also argue that the anisotropic nature of 
steady-state quantities like the shear stress can never- 
theless be captured within such an approximation. 

1. MCT equations for the transient correlators 

The isotropic approximation consists of the following 
three assumptions. First, it is assumed that Fq{t) de- 
pends only on the modulus q — |q|, i.e., 

Fq(t) « F,{t). (Bl) 

Second, we introduce a corresponding approximation for 
the transient cross correlator H^{t) formed with current 
density fluctuations. Since H^{t) is a vector correlator 
whose orientational dependence comes also from the de- 
pendence on A, one cannot introduce such a simple ap- 
proximation like H^{t) « H^it). Instead, we assume 
that the following relation, valid for isotropic quiescent 
systems, to hold: 

The third assumption concerns the modulus of the ad- 
vected wave vector q(t) [see Eq. ([56|l]: 

q{tf =q^ + 2(7f) q^qy + {^tfql (B3) 

We assume that q{tY can be approximated by its ori- 
entational average. This is equivalent to neglecting the 
anisotropic term qxqy and approximating q^ by q^/3 in 
Eq. leading to 

q(i)^q^l + {W/i^m- (B4) 



In the following, we shall see consequences of these as- 
sumptions. 

It follows from the approximation (jB2p 

=^q•Hq(^), (B5) 

implying that q ■ Hq(i) also becomes an isotropic quan- 
tity. Equation (|B5[) is consistent with Eq. (|141ap under 
the isotropic approximation. To see this, we first rewrite 
Eq. (|141ap as 

where the specific form = jSxxSfj.y for the shear-rate 
tensor has been used. The application of the approxima- 
tion (|B1|1 then yields 

|F,W = 7^|F,(0+zq-HqW, (B7) 

because dFq{t)/dqy — {qy/q)dFq{t)/dq. Since now the 
left-hand side depends only on the modulus q, the orien- 
tational averaging of this expression gives Eq. (jB5|) . 
From a partial time derivative of Eq. (jB5[) , we obtain 



g^Fq{t)=^cl■-U^{t). (B8) 
Substituting Eq. (jl41bp into the right-hand side yields 



^Fq{t) = ^igA[q-K.^]iJq\i)-'?'j^qW-E^'?^['^-HqW]'-«[^q-Hq(t)] 



9q 



E/ dstqxM^'^{s)H^,Jt-s)-jY. f dsqxL^^{s)F^^s){t-s). 



(B9) 



Applying the approximation (jBTj) to the second and sixth terms, ()B2|) to the fifth term, and (IB5|) to the fourth term, 
we obtain 



o 2 
A ^ A 

A,M -^0 A 



(BIO) 



where the dot denotes the partial time derivative. The first and third terms on the right-hand side of this equation 
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can be manipulated as 



oqy 



dq q 



(Bll) 



(B12) 



in deriving which we have used the approximations (|B2|) and (jBSp . Both of these terms are anisotropic, and vanish 
after taking the orientational average. [Remember that the left-hand side of Eq. (|B10[) depends only on the modulus 
q.] We therefore obtain 



(B13) 



where we have employed the approximation (jB4p for the modulus of the advected wave number in the fourth and fifth 
terms, and introduced 



F,it) + q' — F,{t) + aF,{t)+ / dsAf-°(s)F^(,)(t-s)+7 / ds V^°{s) Fg^,j{t - s) = 0, 



X.jj, 



(B14) 
(B15) 



Now, we are left with the kernels M^°{t) and iq°(t) which still depend on the wave vector. From Eqs. (jB14p and 
P42ap . one gets for M^°{t) under the approximations (|B1[) and (|B4p 



pv 



2(2^)3 92 [i + (^t)2/3] 



rfk [q • kcfe + q • pcp] [q(t) • k(t)cfe(,) + q(t) • p{t)Cp^t)] Fk{t)Fp{t). (B16) 



It is clear from this expression that the wave-vector dependence of the this memory kernel stems from the terms 
q(t) • k(t) and q{t) ■ p{t). Using Eq. the explicit expression for the former reads 



q(t) • k{t) = q • k + ijt){q^ky + qyk^) -f {jifq^K 



(B17) 



and a similar expression holds for q(t) •p(i). With the same spirit as in the approximation (|B4p . the anisotropic terms 
qxky and qyk^ shall be neglected, and qxkx approximated by q • k/3, i.e.. 



q(t)-k(i) «q.k[l + (7i)V3]. 



(B18) 



Substituting this and a similar approximation for q(t) • p(t) into Eq. (|B16p yields the following expression which now 
depends only on the modulus q: 



pv 



2(27r)3<72 



dk [q • kcfc + q • pc^] [q • \^c^^) + q • pcp(t)\Fk{t)Fp{t). 



Concerning one gets from Eqs. (|B15p and (|142bp under the approximation (|B1 



dk [q ■ kcfc -I- q • pcp] 



2(27r)3 

In this expression, the wave vector dependence comes from 

kxkyit) — kxky -\- {'jt)k^ 



kxkyjt) 'S'fc(f) ^ PxPyjt) ^p(t) 

k{t) Sk(t) p{t) Sp(^t) 



Fk{t)Fpit). 



(B19) 



(B20) 



(B21) 
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and pxPyit). Here again, the anisotropic term k^ky shall be neglected, and fc^ approximated by k"^ /S, i.e., 

k,ky{t) « mk^/3, (B22) 

and PxPyit) shall be approximated similarly. Along with the approximation (|B4p . one then obtains the following 
expression which now depends on the modulus q only: 

„2 A.-, r r S'r, 

Fk{t)Fp{t). 



LTit) = 



it 



2(27r)3 3^1 + (7t)V3 



dk [q • kcfc + q • pcp] 



k „ + p- 



it) 



Pit) 



(B23) 



2. Steady-state quantities 

Under the isotropic approximation (jBip for the transient density correlators and (jB4p for the modulus of the 
advected wave vector, the MCT expressions (|138p and (|139p for the steady-state density correlator and structure 
factor are given by 

F-it) = F,it)+^ r,, [^-%y(^ + ^)^g]^f^-) ^^(, + ,)^^^^^(,)^ (B24) 
"'0 91* + -^gct+s) 







5- ^ — / dq 5- - 5, + / ds --^=AI^== ^ F.isY, (B27) 



From a numerical point of view, it is not necessary to further simplify these expressions since these steady-state 
quantities are the final output of the theory rather than the ones involved in the self-consistent calculations. Of 
course, it is instructive to consider their averages over the orientation q = q/g 

/» 3^/1 + (7»)V3 

to which only those terms in Eqs. (|B24[) and (jB25P proportional to contribute. However, it is more informative 
to regard Eqs. (|B24p and (jB25P as the approximate expressions in which the anisotropic nature of the steady-state 
density fluctuations is retained to the lowest order: such anisotropy arises from the terms in Eqs. (|B24p and (jB25P 
proportional to qxqy Indeed, such a viewpoint is necessary to correctly understand "the isotropic approximation for 
the steady-state shear stress" , adopted in Ref. [l^ , which sounds contradictory since the shear stress is intrinsically an 
anisotropic quantity and vanishes under isotropic density fluctuations. We shall come back to this point in a moment. 

Under the isotropic approximations (jBip and (|B4p . one obtains from the MCT expression (|134p for the steady-state 
shear stress Css 



jW r , /• klky{ky+^Skx)S'kS'^ 

2(2^)3 X ''V kk{s) SIS^ 



ds I dk '^-'^'^V^;;"'-" J,;,'" F.isf. (B28) 



One easily understands that only the term proportional to k^ky survives after the integration over the orientation 
k = k/fc, yielding the following expression for CTss under the isotropic approximation: 

Tfrmm Jo ' ^ "^'^'^ ■ ^^''^ 

This is essentially the same expression as adopted in Ref. [l^ . 

To connect such an isotropic expression for ergs with anisotropic density fluctuations, we rewrite the term in Eq. (jB28|l 
which survives after the integration over k in the following form 

ksT f kxky Sf^ f ■ f , kxky ^k{s) 



The quantity in the curly brackets is exactly the afore- mentioned anisotropic term in Eq. (|B25p . Thus, the 
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steady-state shear stress riss can be handled within the k^ky and the other from the anisotropic part of the dcn- 
isotropic approximation since its MCT expression takes sity fluctuations, which altogether behaves as an isotropic 
a form of the product of two anisotropic terms, one from term inside the integral. 
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